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Abstract. We present a graphical interface designed to demonstrate the techniques 
of radio interferometry used by telescopes like ALMA, e-Merlin, the JVLA and SKA, 
in a manner accessible to the general public. Interferometry is an observational tech- 
nique used by astronomers to combine the signal from a few to tens to hundreds of 
individual small antennas to achieve high resolution images at radio and millimetre 
wavelengths. This graphical interface demonstrates how the number of antenna, their 
position relative to one another and the rotation of the Earth allow astronomers to 
create highly detailed images at long wavelengths. 

-Accepted for publication in the European Journal of Physic^ 



1. Introduction 

The technique of interferometry is widely used in the field of radio astronomy, with 
some of the worlds largest scientific endeavours based upon it (e.g. the Atacama Large 
Millimeter/ Sub-mm Array[H[2] and the up coming Square Kilometre Array[3j Hj). Over 
the coming years telescopes such as these will revolutionise astronomy and hopefully 
inspire a new generation of students to take an interest in science. With this in mind 
the need to explain the techniques of interferometry beyond the most basic level becomes 
apparent, however a complete treatment of these techniques requires an understanding 
of several areas of mathematics beyond that taught prior to undergraduate level, 
namely Fourier transforms and complex numbers. For completeness we have included 
in this paper such mathematics and its application to interferometry, we aim to have 
presented this in such a way that the paper is accessible to readers part way through 
an undergraduate Physics degree (e.g. 2 nd year and higher) and beyond. 

| This is an author-created, un-copyedited version of an article accepted for publication in European 
Journal of Physics. IOP Publishing Ltd is not responsible for any errors or omissions in this version 
of the manuscript or any version derived from it. The Version of Record is available online at 
doi:10.1088/0143-0807/34/l/7 
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Whilst these areas of mathematics are vital for an interferometer they are required 
only in the 'intermediate' stage of interferometric data processing, the stage in which 
the signals received at each individual antenna are combined to produce the final image. 

To this end we have developed a graphical tool to demonstrate some of the real 
life effects on an interferometric image of antenna position and number and Earth 
Rotation synthesis, without outwardly discussing the mathematics required at the 
intermediate stage to generate the final image. The tool itself is written in the Python[§] 
scripting language and can be run on Windows, Mac and Linux. The tool is called the 
"Pynterferometer" , being that it is a Python based interferometer. 

This paper is organised as follows; first we address why astronomers use 
interferometry, and then present the fundamental equations used in the derivation of a 
final image from an observation and the effect of increasing the number of antennas. The 
layout of an array and the use of Earth's rotation in creating a high quality observation 
are then discussed. The approach used by our graphical tool is then discussed and 
an outline of the user interface and example demonstrations are outlined. Finally the 
reception of the tool by the general public is discussed. 



2. Why use interferometry? 

For a conventional telescope the maximum resolution is approximately described by: 

= 1.22^ (1) 

where 9 is the angular resolution in radian^jjj A is the wavelength and D the diameter 
of the mirror/parabolic dish. For optical wavelengths (A ~lCT 7 m) the mirrors required 
to produce angular resolutions of milliarcseconds need only be a few tens of meters 
in diameter. However, in the mm/radio regime, with wavelengths A ~10~ 4 - 10 _3 m, 
antennas would need parabolic dishes of ~ 10 4 m to 10 7 m. The largest steerable single 
dish radio telescope have diameters of up to 100m, beyond this major feats of engineering 
(not to mention finance) are required. 

So in order to achieve resolutions comparable with optical telescopes the field of 
mm/radio astronomy instead uses the technique of interferometry, where the signals 
received by two or more antennas are combined together to approximate having a single 
large dish. The maximum resolution of an interferometer is described by: 

* « ~ (2) 

-E>rnax 

where the symbols have the same meaning as Equation [I] but now D has been replaced by 
B max the maximum separation, or baseline, between two antennas of the interferometer. 
The larger the longest baseline the greater the achievable resolution. 

§ see python . org 

|| 9 is typically referred to in arcseconds and can be converted from radians by multiplication by 

180x3600 



A graphical tool for demonstrating the techniques of radio inter ferometry 



3 



3. A basic interferometer 

The simplest case to consider when describing an interferometer is one consisting of 
only two antennas. In this section we follow the analysis of such a system by [5] and 
[0] . Figure [I] shows the geometric setup of a two antenna interferometer and the rays of 
an incoming astronomical signal. 




Figure 1. Diagram of basic two element interferometer. Blue lines represent the 
incoming wavefront of photons at times t and t + T geo . 



Consider photons arriving in a plane wavefront at our two antennas from an 
astronomical object, with sky brightness distribution J; the wavefront will be incident 
on antenna x before antenna y. With the time delay between the wavefront arrival at 
x and y given by: 

b • s bs cos p . . 

Tqeo = = 

C C 

where b is the baseline length in metres, s a unit vector in the direction of the 
astronomical source, p the angle between b and s, and c the speed of light. This delay 
is known as the geometric time delay as its existence is entirely due to the geometry 
of the system and the finite speed of light. The signals out of each antenna are of the 
form: 

x(t) = V\ cos27rz/t (4) 

y(t) = v 2 cos 2nu(t + r geo ) (5) 

where i>i,t>2 are the signal voltage from each antenna, v the observing frequency and t 
time. The signals from each dish are fed into a computer known as a correlator where 
the signals are cross correlated, (denoted by ®): 

R X ,y{Tgeo) =X®y = x(t) * y(-t) = X(u)Y*(u) (6) 
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here * denotes the convolution of signals x(t) and y(t) (with a time reversal for signal 
y(t)). This convolution is identical to the multiplication of the Fourier transform of x(t) 
denoted X(u) with the complex conjugate of the transform of y(t) denoted Y*{y). This 
ultimately gives us a correlator output: 

Rx,y{Tgeo) = V]_V 2 COS 27lUT geo (7) 

As we are observing astronomical signals the voltage outputs of antennas x and y, 
(^1,^2), are related directly to the brightness of the astronomical object, /(s), over a 
solid angle dfl "seen" by our telescope and the area of the each dish we use to observe 
it, A(s). Leading Equation [7] to become: 

Rx, y ( T geo) = Av J A(s)I(s) cos 2nh\ ■ s dVl (8) 

with b\ being the baseline length measured in wavelengths of the observations (i.e b in 
metres divided by A). Radio/mm telescopes do not observe at a single unique frequency 
but instead cover a number of adjacent frequencies, known as a bandwidth, this is 
symbolised as Av. 

As an observation made with an array telescope contains information from not just 
a single point but an area of the sky (the solid angle denoted by dfl), a reference point 
on the sky known as the phase tracking centre is used within the observations. This is 
the point along vector So, with all other observed points at offset vectors a such that 
s=s +<r. For the phase tracking centre an instrumental time delay is introduced to 
make r geo = meaning all time delays in the signals are due to the positional offset a. 

Finally we define a function known as the complex visibility as, 

V = \V\e^ v = J A(a)I(a)e- i2whx - a d£l (9) 

where A(a) indicates we have normalised, by a factor A , the antenna responses within 
the array and all other symbols have the meaning previously described. Substituting 
Equation [9] and the vector addition, s=s +o", into Equation [8] we are left with, 

R x , y = A \V\ Au cos(27rb A ■ s - (j> v ) (10) 

representing the correlator output. From which we can recover the unknowns which are 
directly related to the source, \V\ and 4>Vi by Fourier transform of this value to obtain 
a component of the sky brightness distribution /(cr). 

3.1. A convenient coordinate frame 

For interferometric observations there is a convenient coordinate frame within which to 
calculate the complex visibility, which is key to the function of the Pynterferometer. 
We define the coordinate w along vector s, with a plane perpendicular to w defining 
the east- west coordinate, u, and the north-south coordinate, v. Figure [2] represents this 
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graphically. As, in this frame, the offset vector a is parallel to the uf-plane it can be 
defined in u,v coordinates, giving Equation [9] as, 

V = I ' A{l,m)I{l,m)e- l2 ^ ul+vm ^ n dl *™ 2 (11) 



m 



where / and m are directional cosines of the vector s. These can typically be replaced 
with the actual angular offsets (i and j) using the small angle approximation (owing to 
the small angles subtended by I and m). So V(u, v) is then simply the Fourier transform 
of the sky brightness distribution of the target observation. 




3.2. Number and positions of antennas 

From the above mathematical approach we can see that each pair of antennas for each 
time interval of observation gives a complex visibility measurement at a single u, v point 
(plus its complex conjugate point). To recover the optimal amount of information about 
an astronomical object observation at as many u, v points as possible are desirable. For 
an interferometer of iV antennas there are iV(iV — l)/2 antenna pairs and therefore 
N(N — 1) u, v points in a single measurement. With each antenna pair acquiring 
spatial information of the scale set by the baseline separation of the antennas (following 
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Equation [2J). So the greater the number of antennas in an array, the greater amount of 
information can be recovered about the target. Plus the greater the number of unique 
baseline lengths the larger the number of spatial scales your interferometer can sample 
in the target object. 

For interferometers built with a small number of antennas, e.g. the Australia 
Telescope Compact Array (ATCA) [7J with 6 antennas or e-MERLIN [8] with 7 antennas, 
the number of u, v points attainable in an instantaneous observation is very small. To 
overcome this a technique known as Earth Rotation Synthesis can be employed. As 
the Earth rotates the antennas in an array will effectively rotate around the vector s 
observing new l,m points ergo a new V(u,v) measurement. By making multiple short 
observations of the same object at different times throughout a days observation, the 
amount of unique u, v points can be increased thus gaining more information on the 
target object. 

To illustrate this each panel in Figure [3] shows the affect of increasing the number 
of antennas in an array (2, 5 , 10, and 50 antennas, top to bottom) and in increasing 
the amount of time on source (instantenous image, 2, 4 and 6 hours, left to right). We 
can see from this that short observations with a large number of antennas over a large 
number of different baselines will achieve a greater total number of unique u, v points 
than several hours observation time with a small number of antennas. This is the reason 
modern interferometers (such as ALMA, MeerKat, ASKAP and ultimately the SKA) 
all make use of large numbers of antennas. 

*** 

Note: In this section, for mathematical simplicity we have avoided several steps used in 
the conversion of antenna signals to a sky brightness map. For an in depth and rigorous 
description of interferometry we again refer the reader to [5] and [B] which were used 
extensively in the writing of this section. We also refer the reader to two key textbooks 
in the subject of interferometry and radio astronomy which are [0] and [TU], both of 
which present a detailed discussion of interferometry. 

4. Simulator approach 

As the primary goal of the Pynterferometer is to demonstrate the affect on observation 
by changing the position and number of antennas within an array, a large portion of the 
mathematical and real world intricacies of an actual interferometer are avoided. The 
Pynterferometer takes two inputs; the user defined array configuration and an input 
image (in .png format) and will give one or two outputs depending on the version 
used. The 'Basic' version outputs the simulated output image, whereas the 'Advanced' 
outputs both the simulated output image and a plot of the wt>-coverage created by the 
user defined array. 

Given a telescope array defined by the user, the Pynterferometer calculates the u 
and v points this array would observe and outputs these as a unit matrix. Next, the 
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Figure 3. it, ^-coverage of an interferometer set out in a logarithmic spiral pattern 
comprised of 2, 5, 10 and 5 antennas (top to bottom) and observing for 10 seconds, 2, 
4, and 6 hours (left to right). 

input image is Fourier transformed using a Fast Fourier Transform algorithm and the 
Fourier components matrix multiplied by the observed uv matrix. This sets all the 
points in the transformed input image without a corresponding uv point to zero but 
retains all information where there is a point. This final matrix is then inverse Fourier 
transformed back into the spatial domain and converted in to an image format resulting 
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Figure 4. Schematic of the dataflow and operations conducted throughout the 
simulation. 



We make a few other assumptions with our approach that a radio astronomer would 
not be able to make. These are that no instrumental or atmospheric noise is included 
and no deconvolution of the image (as is routinely done in practice such as CLEANing 
[TTj ) are undertaken. 

5. The Pynterferometer 

5.1. The interface 

Figure [5] shows a screen grab of the basic version of the Pynterferometer. The gridded 
area (top left) with blue circles is the "Array Configuration" area, users can drag the 
blue circles which represent antennae in the the array to any position within the grid. 
The array configuration in this area will then be used to calculate the output image. 
The input image is to the right of the 'Array Configuration" area and is the image which 
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is used as the "sky" image for transform. The output image (top right) is the resulting 
image of the transform using the input image and the defined array configuration. 
Example output images can be seen in Figure |6j each image was made with a different 
antenna configuration. 

The lower portion of the Pynterferometer consists of a number of buttons, these 
(in columns from left to right) are 1) preset antenna configurations, 2) preset antenna 
numbers, 3) antenna addition and removal buttons, 4) image selection, either to use a 
webcam or one of the two image presets, 5) observing mode selection and 6) antenna 
area selection. Finally on the right are 3 status panels which describe what the 
Pynterferometer is doing (i.e. calculating uv coverage) and the current number of 
antenna and the original configuration of the array. Two final buttons exist below 
the status boxes "Scientific" and "Colourful" which change the colour scale used in the 
output image. The bottom right of the Pynterferometer then includes the logos of the 
home institutions of the authors. 
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Figure 5. A screen grab of the Pynterferometer running on the Mac OSX operating 
system. The images and buttons are described in the text. 
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Figure 6. Example Pynterferometer output images using the webcam input image 
mode and "Scientific" and "Colourful" colour scale settings. Each of these six images 
was made with a different antenna configuration, the effect of which can be seen image 
to image. 

5.2. System Requirements and Availability 

The Pynterferometer works on Windows 7 and Vista, Mac OSX and Linux operating 
systems provided the appropriate Python modules and libraries are provided. The 
prerequisites are packaged up with the Pynterferometer and can be downloaded from 
http : //www . jb . man . ac . uk/ pynterferometer/. 

6. Example demonstrations 

The Pynterferometer can be used to quickly demonstrate the effects of antenna 
configuration/numbers/separation on the output image. Whilst we provide these two 
examples, there are many ways in which the public or a demonstrator can interact with 
this software to learn/teach about the imaging capability of interferometric telescopes. 

6.1. Example 1: Number of antennas 

Using the Add/Remove antenna buttons you can show how increasing the number of 
antennas improves the output image. 

(i) First click 'Spiral' and 'Snapshot' 
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(ii) Click 'Single Dish'. The output image will become a blob (a 2D Gaussian convolved 
with the input image in reality). 

(iii) Click 'Add antenna', the output image will then become a series of stripes. Next, 
increase the number of antennas using the 'Add antenna' button or the presets 
'05/10/30 Antennas', the image will slowly improve in quality. 

6.2. Example 2: Array Separation 

Using the 'Decrease Array Size' , 'Normal Array Size' and 'Increase Array Size' buttons 
you can show how a larger antenna separation will observe more of the small scale 
structure, whilst a closely packed configuration will see only the extended structure. 

(i) First click 'Use Galaxy' and 'Earth Rotation'. 

(ii) Click 'Decrease Array Size', the grid lines in the array panel will increase in 
separation, these grid lines always show boxes of 500m x 500m across so you 
are looking at a small area of landj] . The output image will show the galaxy fairly 
well but most of what is seen is the extended/large scale structure (e.g. the overall 
shape of the galaxy). 

(iii) Click 'Normal Array Size' the grid lines will become more closely packed so you 
are looking at more land (i.e. the antennas are now further apart from one another 
with max separation = 7km). An increase in the amount of smaller scale structure 
seen in the output image will now be seen. 

(iv) Finally, click 'Increase Array Size' again the grid lines will become more closely 
packed so you are looking at yet more land (i.e. the antennas are now even further 
apart from one another; max separation = 14km). Here the output image will be 
devoid of most extended structure and dominated by fine small scale structure. 

7. Reception of the tool by the public 

This tool was originally created for use on the ALMA display at the Royal Society 
Summer Science Exhibition 2012, which ran from the 3 rd to the 6 th of July 2012. The 
event saw 11,090 visitors during the 6 days it was open with a large proportion passing 
and speaking to the staff on the ALMA stand. Not all stand visitors used the tool but 
we estimate over a thousand uses over the whole event. 

During the exhibition there was a twitpicj^] feature, which posted the users image to 
the URL http://twitpic.com/photos/alma_sse/ was included. Over the exhibition 
week the 'twitpic' feature, saw ~200 uses. Please note this feature is not available in 
the main Pynterferometer release. 

Some members of the public were happy to use the tool unsupervised and asked 
questions when they felt it necessary, whilst the majority were lead through the tool by 

% Please note the telescope circles are oversized for ease of use. 
+ twitpic.com 
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one of the stand staff. The staff found this as a very useful tool for visually demonstrating 
the process of radio interferometry. It was used quite heavily when trying to overcome 
the difficulties of exampling why the Earth rotation helps. 

8. Summary 

We have presented a graphical interface designed to demonstrated the techniques of 
radio interferometry in a manner accessible to the general public. The Pynterferometer 
demonstrates how the number of antenna, their position relative to one another and 
the rotation of the Earth allow astronomers to create highly detailed images at long 
wavelengths. The aim of the Pynterferometer is to introduce the observational technique 
of interferometry which will become increasingly prevalent in professional astronomy in 
the next few years as e.g. ALMA, ASKAP, MeerKAT and the SKA all become fully 
functioning instruments. 

Acknowledgements 

The authors would like to thank Dr John Richer who lead the ALMA team at the Royal 
Society Summer Science Exhibition 2012 and thought this would be a good idea! The whole 
team who manned the stand at the Exhibition, the staff of the UK ARC Node and PDRAs 
and PhD students from both Manchester and Cambridge who tested this software. The Royal 
Society for selecting the ALMA team to display at the Exhibition. The STFC as well as 
the University of Cambridge, The Kavli Institute for Cosmology and the UK ARC Node for 
funding the ALMA Exhibition. 

[1] ALMA. ALMA observatory public website, "http://www.almaobservatory.org", August 2012. 
[2] R. L. Brown, W. Wild, and C. Cunningham. ALMA - the Atacama large millimeter array. 

Advances in Space Research, 34:555-559, January 2004. 
[3] Square Kilometre Array. SKA telescope website, "http://www.skatelescope.org", August 2012. 
[4] R. T. Schilizzi, P. E. F. Dewdney, and T. J. W. Lazio. The square kilometre array. In Society of 

Photo-Optical Instrumentation Engineers (SPIE) Conference Series, volume 7733 of Society of 

Photo-Optical Instrumentation Engineers (SPIE) Conference Series, July 2010. 
[5] G. B. Taylor, C. L. Carilli, and R. A. Perley, editors. Synthesis Imaging in Radio Astronomy II, 

volume 180 of Astronomical Society of the Pacific Conference Series, 1999. 
[6] Bernard F Burke and Francis Graham-Smith. An Introduction to Radio Astronomy. Cambridge, 

2nd edition, 2002. 

[7] Australia Telescope National Facility. Australia Telescope Compact Array website. 

"http://www.narrabri.atnf.csiro.au/", September 2012. 
[8] e-MERLIN/VLBI National Facility. e-MERLIN/VLBI National Facility website, "http://www.e- 

mcrlin.ac.uk/", September 2012. 
[9] A. R. Thompson, J. M. Moran, and G. W. Swenson. Interferometry and Synthesis in Radio 

Astronomy. John Wiley & Sons., 2nd edition, 2001. 
[10] K. Rohlfs and T. L. Wilson. Tools of Radio Astronomy. Springer., 4th rev. and enl. ed. edition, 

2004. 

[11] T. Cornwell, R. Braun, and D. S. Briggs. Deconvolution. In G. B. Taylor, C. L. Carilli, and R. A. 



A graphical tool for demonstrating the techniques of radio inter ferometry 13 

Perley, editors, Synthesis Imaging in Radio Astronomy II, volume 180 of Astronomical Society 
of the Pacific Conference Series, page 151, 1999. 



